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We investigate collision dynamics involving two electrons and a quantum dot embedded in a quantum wire. 
One electron is initially at rest in the ground state of the dot, whereas the other electron is incident on the dot 
with a rather well denned velocity. Excitation, reflection and transmission probabilities are obtained by resolving 
the dynamics explicitly. Also, the probability of both electrons being unbound after collision is studied. It is 
found that the dynamics is strongly influenced by the presence of doubly excited states; reflection becomes 
considerably more probable in the energetic vicinity of such states. This, in turn, contributes to an increase in 
the probability of trapping also the second electron within the dot. 

PACS numbers: 32.80.Hd, 73.21.-b, 73.63.-b, 33.35.+r 

I. INTRODUCTION 

Since the realization of the first quantum dots about two decades ago, much effort has been spent studying them. This is 
motivated by the fact that such devices hold promise for several useful applications. Also, the possibility to design particular 
systems with specific features makes the study of such structures an attractive one. Although these so-called "artificial atoms" are 
similar to atoms in the sense that they share many features, such as, e.g., energy quantization and shell structure, their structure 
may often be tuned while the structure of atoms is fixed by fundamental constants. This opens up a wide area for studying 
quantum phenomena. Of course, another feature which is appealing to a theorist, is their reduced dimensionality. 

A number of phenomena well studied in atoms have later been investigated in quantum dots as well. One such example are 
doubly excited states embedded in the continuum of a quantum dot. Such states, which are unbound but in some respects resem- 
ble bound states, are examples of resonances states. Their existence have been widely studied, experimentally and theoretically, 
in all areas where quantum mechanics apply - also in the context of quantum dots, see, e.g., Refs. 1 5 . This is also the aim of the 
present work. We intend to do so by resolving the explicit time evolution of the system. Although resonances are usually studied 
in a time independent framework, other dynamical examples are also found in literature, see, e.g. RefsP^. 

Our system consists of two electron in which one is free and the other is trapped in a quantum dot initially. As in, e.g., RefsPEH 
we will take our system to be confined in two spatial directions so that our problem becomes an effectively one-dimensional 
one, i.e., the quantum dot is embedded in a quantum wire. In a collision event like this, there are several possible outcomes. 
As the system is unbound, at least one electron will escape unless the system is allowed to lose energy somehow. The escaping 
electron(s) may either be reflected or transmitted. If one electron remains bound, it may be found in an excited state, resulting in 
an inelastic scattering event. If the the projectile is incident with sufficient energy, it may liberate also the bound target electron, 
analogous to ionization in atomic and molecular physics. Here we will refer to this process as knockout. 

One dimensional scattering problems have been studied extensively - also with particular attention paid to the role of res- 
onances, see, e.g., Refs! 13 l 14 l Such studies tend to involve only one particle, i.e., the resonances are shape resonances. In the 
present work we will focus upon how the doubly excited states which in fact are present in our system influence the reflection and 
transmission probabilities. We will also investigate the capture dynamics in this regard. In Ref. 15 , which also concerns electron 
capture and doubly excited states, we may read the following: "An interesting question is whether and how a second electron 
can be captured when a first one is already present in the dot." This is precisely what we address. In RefP^two processes are 
studied: single electron capture by emission of a phonon, and relaxation of one electron via intraband Auger decay, i.e., two 
electrons in a doubly excited state decay into a state in which one electron is relaxed and the other one escapes. In this process, 
no phonon is emitted; the excess energy is carried away by the Auger electron. In the present work we aim to take advantage 
of both the (intraband) Auger effect and spontaneous emission in order to capture the projectile electron in the same dot as the 
target electron. In fact, this process is the completely analogous to what is called dielectronic recombination in atomic physics, 
which, in turn, is the time reverse of Auger decay. 

Although our starting point is the Schrodinger equation, the equation we actually solve in the end is the Lindblad equatiorl^-^. 
Usually, master equations are used to accommodate for irreversible processes such as, e.g., dissipation. Moreover, in the con- 
text of quantum transport, they allow for a convenient description of the coupling between nanostructures and their leads 18 19 . 
However, this is not the motivation in the present work. It is simply a consequence of our desire to impose absorbing boundary 
conditions. 

We have chosen to work with units defined by letting h, the effective electron mass m* , the elementary charge e and 
l/(47T£ r £o), with e r and e being the relative permittivity and the permittivity of vacuum, respectively, define the units of 
their respective quantities. For, e.g., GaAs this corresponds to a length unit of about 10 nm, a time unit of about 60 fs and an 
energy unit of 1 1 meV. These units will be used throughout. 
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The paper is organized as follows: In Sec. [II] the theoretical framework applied is presented and, in part, derived. The results 
are presented and discussed in Sec.|IIIJ while conclusions are drawn in Sec.|IV| 

II. THEORETICAL FRAMWORK 

The dynamical equations which are to be solved in this work, will deviate considerably from our starting point, which is the 
Schrodnger equation for two interacting particles in three dimensions. After having reduced the problem to an effectively one- 
dimensional one, we will, for numerical convenience, introduce absorbing boundary conditions. This calls for the introduction 
of the Lindblad equation. Finally, we demonstrate how we estimate the probability of capture through spontaneous emission of 
a photon or a phonon via the Fermi golden rule. 



A. Effective potentials 



The dynamics is governed by the time dependent Schrodinger equation, which in our units reads 
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where V is some local confining potential, i.e. the quantum dot, and W is the electron-electron interaction. 

We will restrict the dynamics to one spatial dimension, i.e., the system is confined to the ground states in two Cartesian 
directions, which we will take to be the y and the z direction. We will take the confinement to be represented by infinite square 
wells of equal extension, /, in both these directions. Assuming that the confinement is strong enough for us to neglect correlation 
in these directions, we may, up to a trivial phase factor, write the wave function on the form 1 ^! 



$(rx,r 2 ;t) = ^(x 1 ,x 2 ;t)ip (yi)ipo(zi)ipo(y2)ipo(z2) with 
ip (x) = Jycos(yi), \x\<l/2 . 
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The assumption that the system remains in the ground state in the y and z directions requires that the typical energy transfers 
involved in the processes we wish to study, AE, is considerably smaller than the energy it takes to excite the system in the y or 
z direction, i.e., AE < 3n 2 /(2l 2 ). 

The problem is reduced to an effectively one dimensional one by integrating out the y and z dependencies: 
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The effective interaction may be written aP 
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which may be effectively calculated numerically by means of, e.g., Monte Carlo integration. This interaction potential is plotted 
in Fig.[T]along with the potential 

W(xi 2 ) = —=2=== . (11) 

\/Xl 2 + (IS) 2 
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In this work, as in several others such as, e.g., RefsP^l the latter has been applied in the calculations. This is motivated by the 
fact that this potential has a continuous derivative, whereas W features a cusp at x\2 = 0, c.f. Ref. 11 . The asymptotic behavior 
is the same for both potentials, and they scale with the confinement I in the same manner 9 . The parameter S is chosen such that 
the difference between W and W integrates to zero. This is achieved with S « 0.275. We do not see any a priori reason why 
replacing W with W should alter the phenomenological features of the system. 

As in Ref. 3 , the effective one-particle potential, V, is chosen to be a negative Gaussian: 

V(x) = -D v e X p(-^j , (12) 

where Dy is the depth and ay provides the width of the potential. From a theoretical point of view, this is a convenient choice 
as it allows for analytical continuation into the complex plane, which will be addressed shortly. Moreover, the short range nature 
of the potential ensures that there are no Rydberg states. 



B. The spectrum 

As pointed out initially, the importance of doubly excited states is crucial to this work. Thus, we need to be able to localize 
these resonances. This is readily done by the standard technique of complex rotatiorPSEl our spatial variables are multiplied 
by a complex phase factor, x — > e l0 x, with < 9 < tt/4. The resulting time independent Schrodinger equation is solved 
while maintaining Dirichlet boundary conditions. As is well known, this provides "eigen energies" which are unaltered for 
bound states and rotated the angle 29 into the fourth quadrant of the complex plane for "ordinary" continuum states. The doubly 
excited states, however, are identified by the fact that their complex eigen values are virtually independent of 9 (for large enough 
9). The real part of such an "eigen energy" gives the position of the resonance, whereas the imaginary part provides its width. 

An example of such a complex spectrum is shown in Fig. [2] It corresponds to the system at hand with the parameters Dy = 4, 
ay = 1.5 and I = 0.7. The eigen energies of the corresponding one-particle system, i.e., the thresholds, are also displayed. 
From this figure we see that the two-particle system features two bound states; the excited state is just barely bound. Morover, 
one eigen energy with non-vanishing imaginary part is seen to be ^-independent. This complex eigen energy is located at 

E tes = -1M7- 0.089m = E Ies -iT Ies /2 (E Ies , T ies G R) . (13) 



C. Resolving the dynamics 



Our initial state is such that one electron is free and another one is trapped in the (one-particle) ground state for the confining 
potential V. The free electron is given by a Gaussian wave packet separated from the bound electron. This is illustrated in Fig. [3j 
The packet is centered around a mean velocity fc with a rather narrow width in momentum space. More specifically, we take 
the spatial part of our initial state to be 
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xq is the initial centering of the incoming electron, and to is the instant a freely propagating wave packet is at its narrowest 
(spatially). 4>gs(x) is the one-particle ground state. As seen from the exchange symmetry of the above state, we have chosen our 
initial state to be a spin singlet state. The mean energy of the incoming particle is given by 
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with standard deviation provided by 



The total energy of the system is 
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where the last term is the one-particle ground state energy. 

Since the incoming Gaussian is rather narrow in momentum space, it must be quite wide in position space. Thus, a numerical 
grid of rather large extension is needed. During propagation, this situation is complicated further due to dispersion and the fact 
that parts of the wave packet will be reflected. In terms of feasibility, this may be problematic for any system involving more 
than one particle - even in the one dimensional case. 

Thus it is tempting to introduce absorbing boundary conditions as a numerical tool. However, doing so introduces other 
problems: When absorption takes place, the entire wave function, which is normalized to the probability of having two particles 
on the grid, vanishes. If we want to know what happened to the possibly remaining electron, we are at a loss. We will, e.g., not 
be able to know whether one or two electrons have escaped, or what the excitation probabilities are. As we also aim to calculate 
the probabilities of such events, we are indeed faced with this problem. 

As it turns out, the solution to this problem is provided by the Lindblad equatiorP^l. This is explained in detail in Ref.^1 
Here we will very briefly outline the main ideas. The Lindblad equation is the natural starting point as it ensures conservation of 
positivity and trace for a Markovian process, such as absorption. By comparing the effective Hamiltonian involving a complex 
absorbing potential with the generic structure of the Lindblad equation, a source term is identified, which "feeds" a one-particle 
density matrix pi as the two particle state vanishes due to absorption. The one-particle system, in turn, is also exposed to the 
absorber, which may lead to population of the vacuum-state po = po\— }{— 1> i.e., the state in which there are no particles. The 
total density matrix features "super selection", i.e., if the initial state corresponds to a well defined number of particles, only 
"diagonal blocks" are populated; 
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P = ^2Pn,n , (19) 

ri=0 

where n corresponds to the number of particles. 

In this work, we will modify the theory of Ref. 24 slightly in the sense that we will discriminate between the left and right 
absorber. We will distinguish the one-particle density matrix in a part arising from absorption at the left edge and one arising 
from absorption at the right edge, 

Pl = p\ + pf . (20) 

Similarely, we may distinguish the population of the vacuum state, p , in four parts corresponding to the four ways of absorbing 
both electrons. This is illustrated in Fig.|4] Specifically, in our grid representation, the dynamical equation may be written as 

ip = [H , p] — iCl [p] — iCr [p] with the Lindbladians (21) 



£l/r = 22^h/R(xn) ({4c„,p} - 2c 



nPcl) 



where Cn annihilates (creates) a particle in position x n and r L / R is the left/right absorber. H is the total Hamiltonian expressed 
in terms of second quantization. 

In this work we have chosen to use an absorber of the form 

T R (x) = \ 2 - 5fc ° {X ^Lt^ ' L / 2 ~ ^ abs <X<L / 2 ( (22) 
1 0, otherwise 

where Z a b s = L/10, i.e., the absorber is 'turned on' at distances one tenth of the extension of the grid from the edges, and it 
increases quadratically towards 2.5/cq « 5E at the edge right edge at a; = L/2. The left absorber is defined correspondingly, 
and their sum constitutes the total absorber, r = IY + Tr. We have chosen the strength of the absorber proportional to the mean 
energy of the incoming electron in order to achieve sufficient absorption while inducing minimal unphysical reflection. 

If we write down the evolution of the various constituents of the total density matrix separately, it is seen that the two-particle 
part remains a pure state, p2 = \^ 2 ) (^2 1 where the wave function obeys the ordinary non-Hermitian Schrodinger equation, 

i* 2 = \hi+h 2 + W(x 12 )-iT(x- L )-iT(x 2 )\^ 2 . (23) 

The one-particle parts are given by 

ip? = [h,p?]-i{r,p?} + 2iJ2rD(x n )c n \y 2 )(V 2 \c$ l (24) 

n 

and the vacuum probability is given by 

^ B = 2(-|^r B (x„)c„pfcJ l |-> . (25) 

n 

The indices '_D and stands for either 'left' or 'right' (L/R). 
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D. Decay by spontaneous emission 

This work aims to investigate whether doubly excited states may be used actively in trapping electrons. Of course, as long as 
the dynamics does not involve any interaction with any field, energy is conserved and capture is prohibited. Thus, the possibility 
of relaxation via emission of some field quanta must be included in our calculations somehow. Of course, this may be achieved 
by actually including some quantized field in our formalism. This is numerically infeasible as this would introduce too many 
additional degrees of freedom in our dynamical calculations. Another way would be to describe the field classically, i.e., non- 
quantized, which is admissible if a large number of field quanta are present initially. However, we are interested in spontaneous 
emission in this context. 

With our dynamics already described by an equation of Lindblad form, an approach in which the degrees of freedom of the 
quantized field is traced out, resulting in a Master equation for the reduced density matrix, could be a fruitful pai n 25 ! 26 ! However, 
since our two-particle wave function has a rather well defined energy, we may resort to a simpler strategy: We will estimate 
dynamical capture rates based on the Fermi golden rule. In the limit that the electronic system couples weakly to the field, the 
rate at which the system is captured into bound state b by spontaneous emission of a field quantum is estimated by 



n(t)-27r^|(<f b ,l k |i/ 7 |vl/ 2 (t),0)| 2 6 k 2 j2 + e« - (eg + Hu(k)) 



(26) 



where the sum runs over all modes k. |$t>, lk) is the composite state in which the electronic part is in bound state b and 
one quantum has been emitted in mode k, \^>2(t), 0) corresponds to the electronic state ^ with zero field quanta. Thus, by 
calculating the time integral of Eq. ( |26| , the probability of capture through spontaneous phonon emission may be estimated. 
We should also address the interaction Hj. The interaction with the phonon field may be represented by 



photon = ir . £ ^ ( ek>k ^x.r „ e^a^e-*") , (27) 

where V is the quantization volume and annihilates (creates) a photon in the mode given by the wave vector k and polar- 
ization fi. However, in a solid state context, we would usually expect the the interaction with phonons to be more important. In 
the case of the deformation potential, this interaction reads 

rrphonon \ \ "> k / - ikr . ~t — ikr\ /"10\ 
H I =A Z^^T7Tl ake +a k e j ' (28) 

k V w O) V ' 

where the parameter A depends on the volume of the sample, the density of the material and the strength of the deformation 
potential tensor 2 ^!. 

Although seemingly rather similar in form, these interactions have quite different features. In the case of the photon interac- 
tion, the wave number in the exponential is related to the energy shift to by k — lu/c with c being the speed of light. As this 



is a large number and u> is typically moderate, we may safely approximate the exponential factors in Eq. ( 27 1 by unity, thus 



removing the explicit energy dependence of the interaction. This approximation, i.e., the dipole approximation, is one that only 



couples states of opposite parity. When it comes to the phonon interaction, Eq. ( 28 I, the energy shift u and the wave number k 
is related by dispersion relations, which are specific to the material at hand and the the kind of phonons. On our length scale, 
typical fc-values may become quite large, thus resulting in a very strong energy dependence. 



E. Numerical implementation 

The two-particle wave function is calculated using a rather standard spilt operator technique^. We have split the propagator 
according to the number of particles it acts on: 

* 2 (£ + t) = e - iAr/2 exp(-iW(x 12 )T)e- lAT/2 y 2 {t) +0(t 3 ) with (29) 
A = hi + h 2 - i[T(xi) - T(x 2 )] ■ 

We represent the wave function as a symmetric matrix with indices given by the position of each particle. In this way, the two- 
particle operator is implemented as entry-wise multiplication and the one-particle operators are implemented as left and right 
multiplication, respectively, with the state matrix. 

In earlier implementations of the Lindblad equation for an initial two-particle system with absorbing boundaries, the numerical 
scheme has been suffering from not being manifestly trace-conserving 2429 . Thus, as conservation of probability is crucial, it 
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was necessary to resort to a very small time step r. In the present work this problem has been considerably reduced. We use the 
following scheme, based on a three point formula for integration, to solve Eq. (|24|: 



p?(t + 2r) = e~ 2lBT p? [t)e 2lB " T + (30) 



T ( e -2iBr 



S D (t)e 2lB ' T + 4e- lBT S D (t + r)e lB ' T + S D (t + 2r)) + 0(t 5 ) with 
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B = h-iT and 
S D (t) = 2j2^D(x n )c n \^ 2 (t))(^2(t)\ci . 

n 

As before, "D" is either L for 'left' or R for 'right'. 

The four parts of the zero-particle probability, p BE , are simply obtained by integrating Eq. (25 1. Also this may be done by 
means of the three point formula. Examining to what extent the total trace, Tip = l^l 2 + Trp! + p remains unity at all times 
provides a good check for the accuracy of our calculations. 

We will need the (non-complex scaled) wave functions of all bound states in order to be able to calculate capture rates, c.f. 
Eq. ( po}. T hese states are, along with their energies, determined by propagation in imaginary time, i.e., by replacing r by — it 
in Eq. (|29[). With this substitution, an arbitrary initial state with some symmetry property will exponentially converge towards 
the lowest state within that symmetry, and its energy is provided by the change in norm at each time step. Within a specific 
symmetry, excited states are found by repeatedly projecting away already obtained lower lying states. 



F. Analysing the state 

For large times t, the traces of each of the constituent parts of the density matrix are subject to rather straight forward 
interpretation. The trace of p\, e.g., converges towards the probability that one, and only one, electron has been absorbed 
and that it has been absorbed to the left. I.e., it converges towards the reflection probability R. Similarly, Tr pi provides the 
transmission probability T: 

R = lim Trp\{t) 

t—too 

T = lim Trp?(t) . 

t— > oo 

With an unbound initial state and no field interaction the norm of the two-particle wave function will necessarily converge 
towards zero. Moreover, if knockout is energetically prohibited, the reflection and transmission coefficients may be found 



without calculating p\ explicitly. To see this, consider the evolution of the trace as dictated by Eq. (24 1: 



Trpf = 2(* 2 |[r D (x 1 )+r D (x 2 )]|* 2 )-2Tr[rpf] . (31) 

If we know that pi never overlaps with the absorber, the latter term vanishes and we may find T and R simply by integrating the 
overlap between the two-particle wave function and the respective parts of the absorber as functions of time, 

Rn2 (MtWh(xi)+r h (x 3 )]\^a(t))dt , (32) 
Jo 

with T calculated analogously. 

Although this rather intuitive relation may be useful, a lot more information may be obtained in the framework provided by 
the Lindblad equation. We may, e.g., find the probability of exciting the electron that remains after collision and the probability 
that both electrons are unbound after collision. The latter is simply provided by po in the limit t — >• oo, whereas the population 
of bound one-particle state number b is found as 

^ixcitcW = Tr [\<p b )(<p b \pi(t)] = (<Pb\pi(t)\<Pb) ■ (33) 

As we have seen, all of these quantities may also be found differentially in terms of the direction(s) of the escaping electron(s). 
The knockout probability may also be found in a manner which converges faster in time: 

P ko = lim (p (i) + Tr[7Wi(*)]) , (34) 

t— > oo 



where "Pub projects onto the subspace of unbound one-particle states. The latter formula does not, however, provide any 
information about the direction in which the second electron is ejected. 
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As the absorbers remove all information about escaping particles, the system is bound to lose coherence. This loss may be 
quantified, e.g., by the von Neumann entropy, S = — Tr [p log 2 p], or by the purity, £ = Tr p 2 , which is unity for pure states and 



decreasing with increasing "mixedness". Since p is block-diagonal in term of particle number, c.f. Eq. ( 19 1, this quantity may 
be calculated as the sum of the purity of each block: 



c 



1*9 



Tr p\ 



(35) 



III. RESULTS AND DISCUSSION 

In order to illustrate the collision dynamics, a series of "snapshots" of the particle density is displayed in Fig. [5] The full 
curve is the particle density given by the two particle wave function, i.e. \^2(xi,x; t)\ 2 dxi, and the dashed curve is the 
particle density of the one-particle sub-system. The latter is simply the diagonal of p\. In the absorption free region of the grid, 
the sum of these densities coincide with the particle density obtained from a wave function obtained on an untruncated grid. The 
incoming electron has a mean energy E = 1.8 with = 0.1. In these calculations, a grid extending over L = 150 length units 
have been used, featuring 511 grid points. Converged results were obtained with the numerical time step r = 0.05. 

As is seen, the one-particle sub-system is populated as the two-particle wave function is absorbed. As the system is completely 
unbound, ^ 2 eventually vanishes completely, leaving all population in pi (knockout is not energetically accessible in this case). 
This is also illustrated in Fig. |6| which shows the norm of the two-particle wave function as a function of time - along with 
the traces of p\ and p 1 ^. The instances shown in Fig. [5] are indicated at the x-axis. We have also included the evolution of the 
purity, c.f. Eq. ( [35| , of the system. As this quantity converges towards a value below unity, it is clear that coherence is lost in the 
process, and the final state is not a pure state. 

In Fig.|6]we observe that there is a certain probability of reflection in this collision. Although just barely, this may be observed 
in Fig.[5]as well. In this specific case, the reflection probability is 5.3%. In Fig.[5]it is seen that the escaping electron comes 
out in two "lobes" - corresponding to different group velocities. In Fig.|6]this is manifested in a step-like increase in the two 
one-particle traces. This is due to the fact that the incident electron has enough energy to excite the bound one. For the slow 
lobe, this excitation has taken place, whereas the faster lobe corresponds to elastic scattering. 

In Fig. [7] we have shown how the reflection and transmission coefficients depend on the the energy Ei of the two-electron 
system. Here the one-particle ground state energy is e gs = —3.141. We have also displayed the probability for the system to 
end up in the various one-particle bound states. The parameters of the potential V and interaction W are the same ones as those 
of Fig. [2] The incident electron has — 0.06, c.f. Eq. (15). Converged results where obtained using a grid consisting of 1023 
points extending over L — 250 length units. As in Fig. [5 we have used the numerical time step r = 0.05. 

It is seen that for electrons incident with energies below the energy required for excitation, only elastic scattering takes place, 
i.e. the one-particle ground state is the only populated state after collision. When the inelastic channels open, the populations 
of the excited states become appreciable. The dependence on initial energy is not monotonic, however. Their structure seems 
to be related to the onset of each channel, whereas the populations of the various bound states feature more or less plateau-like 
behavior in between. 

The most striking feature of Fig. [7] however, is the pronounced peak in the reflection coefficient R{Ei), or, correspondingly, 
the dip in T, at Ei w E tes . This is a clear manifestation of the doubly excited state seen in Fig. [2] A careful look, however, 
reveals that the maximum in R does not completely coincide with E rcs \ it is shifted slightly to the left. As is indicated in Fig. [8] 
this may be due to the interference between the resonance and the background continuum. Here we have plotted a close up of 
the peak together with a convoluted Fano profile 30 . The top of the peaks coincide for a q-value of —4.5. This parameter and 
the hight of the peak are the only two parameters used in the fitting. The rest is provided by the energy-width of the incident 
electron, o^, and the resonance parameters E Ies and T Tes , c.f. Eqs. ( 15|13 1. Performing a convolution is necessary as the width in 



the energy distribution of the incident electron is comparable to the width of the resonance, r rcs . A second case with a narrower 
energy distribution, = 0.04, is also included. 



In the expression for the initial Gaussian wave packet, Eq. ( 15 1, the parameter to is the time at which the packet, if allowed to 



propagate freely, is at its narrowest. In all the above cases, this has been set to zero, i.e., the projectile wave packet widens as it 
approaches the target. One may ask to what extent the collision dynamics are altered if we chose a different to- We have repeated 
the above calculations with a to which is such that the projectile is at its narrowest at impact. For the quantities presented in 
Figs. [7] and [8] no differences where seen. 

For projectile energies larger than the ionization potential of the target electron, knockout becomes energetically possible. In 
Fig.|9]the probability of this process is shown as a function of the initial energy for the same case as in Fig. [7] A finite probability 
for knockout is seen for energies beyond zero. It does not increase monotonously with energy. It is also seen that the probability 
stays rather low. For this reason a lower numerical time step had to be used than in the case of Fig. [7] The low probability 
of knockout may be understood from the fact that numerical solutions of Newton's second law for the system suggest that the 
process does not happen at all classically. When the process does have a finite probability quantum mechanically, it may be 
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related to the fact that classically, the projectile electron may gain energy on the expense of the target electron, whereas this is 
impossible in a quantum mechanical context as the target electron is in its ground state initially. 

The fact that reflection becomes strongly enhanced due to the presence of doubly excited states suggest that also the capture 
probability may be increased by such resonances. Since a projectile electron which eventually is reflected typically "lingers" 
around the target electron and the confining potential longer than a transmitted one, there should be more time for the free 
electron to "fall into" the confining potential by spontaneous emission of a phonon or a photon. 

Figure 10 shows how the probability of capture due to spontaneous emission of a photon depends on energy of the initial state. 
As with the reflection coefficient, a clear peak is seen in the vicinity of the resonance state. Also in this case, a shift towards 
lower energies as compared to E ICS is seen. This may, once again, be due to the interference, but it could also be due to the 
fact that as the "overlap time" in general decreases as the group velocity of the projectile electron increases. For this reason, the 
capture probability is an over all decreasing function of energy - except in the vicinity of resonances. 

Also for the phonon interaction, the capture probability is a decreasing function of energy. However, here the effect mentioned 
above is completely dominated by another one. Due to the before mentioned strong explicit energy dependence in the interaction 
Eq. ( |28] >, the probability of relaxation via phonon emission falls off quite dramatically. For this reason, no resonance peak is 
observed for phonons for the case at hand. We do, however, still expect to see such peaks at smaller energy scales. 



In Fig. 11 capture probabilities are displayed for a case given by the parameters Dy — ay = 3 and I = 1. The complex 



energy spectrum is shown in the first panel. It is seen that the system supports 5 bound states and two resonance states between 
the first and second threshold. The next panels give capture probabilities for photon and phonon emission, respectively. For 
the former, this have also been shown for positive and negative parity separately. Two peaks are clearly seen - one for each 
resonance. Moreover, it is seen that the resulting population of bound states reflects the parity of each of the two resonances. 
In this case, a resonance peak is seen for the phonon emission as well. Any peak near the resonance with position close to —2 
is hard to see, however, for the same reason as above. The dispersion relation used here corresponds to longitudinal acoustic 
phonons in the [111] -direction for GaAs^. 



IV. CONCLUSIONS 



Collision dynamics between two electrons in a quantum wire where one of them is initially trapped in a quantum dot have 
been investigated. This has been done within a theoretical framework in which the concept of absorbing boundaries have been 
generalized to apply to many-particle systems. We have seen that, as expected, the collision dynamics is strongly influenced by 
the presence of doubly excited states. Specifically, it is seen that the reflection coefficient features a pronounced peak when the 
energy of the incident electron is such that it is resonant with an excitation to such doubly excited states. It was seen that this 
phenomena may lead to a enhanced probability of capture via spontaneous emission. Thus, by tuning the energy of the incident 
electron, doubly excited states may be exploited in order to facilitate the population of quantum dots. 
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FIG. 2: (Color online:) The spectrum of the complex rotated Hamiltonian for three values of 6. The eigen energies of the corresponding 
one-particle system, i.e. the thresholds, are also displayed. 




FIG. 3: The particle density corresponding to our initial state: One electron, represented by a Gaussian, and one electron trapped in a confining 
potential, which is also a Gaussian (dashed curve). 
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FIG. 4: (Color online:) Schematic illustration of the various ways of reducing the number of particles. We start out with two particles, which 
may be turned into a one-particle system by absorption either to the left or the right. The remaining electron may also, in turn, be absorbed, 
thus populating the vacuum state. 




FIG. 5: (Color online:) The particle density at various times during the collision process. The full (red) curve is the particle density originating 
from the two particle wave function, whereas the dashed (blue) curve is the particle density corresponding to the one-particle sub-system. Also 
illustrated are the complex absorbing potential and the confining potential (thin curves). 
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Time 

FIG. 6: (Color online:) These curves refer to the collision process depicted in the "cartoon" in Fig. [5] The full curve shows the norm of the 
two-particle wave function, the dashed one is the trace of pf ', i.e. the probability of right absorption, the dash-dotted one is Tr p\ multiplied by 
10. The purity of total density matrix is shown as a dotted curve. The instants corresponding to the various "snapshots" in Fig.|5]are indicated 
by circles. 




Energy 

FIG. 7: (Color online:) The transmission and reflection coefficients, T and R, respectively, as functions of the energy of the initial two-particle 
wave function, Ei = feo/2 + e gs (full curves). The dash-dotted curves show the probability for the system to end up in the various one-particle 
bound states. The x-marks on the x-axis indicate the onset of each of these channels - except for the first one, which is the one-particle ground 
state energy. The red circle is the real part of the "eigen energy" og the resonance. 
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FIG. 8: (Color online:) A close up on the peak seen in R(Ei) in Fig.[7](full curve) along with a convoluted Fano-profile (dash-dotted curve). 
Also included are the corresponding curves for a case in which the Gaussian wave packet of the projectile electron has a narrower momentum 
distribution (a^ = 0.04). As in Fig. [7] the circle on the a>axis marks the position of the resonance state. 
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FIG. 9: The probability of both electrons being liberated as a function of the energy of the system. 
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FIG. 10: Estimate of how the probability of capture via spontaneous emission of a photon behaves in the vicinity of the doubly excited state. 
The position of the doubly excited state is indicated by a circle. 
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FIG. 11: (Color online:) a) The complex spectrum of the case at hand for three different scaling angles. Two resonance states are identified 
between the first and second threshold, b) Estimates of the probability of capture via photon emission. The partial contributions from bound 
states of positive and negative parity are also displayed, c) Estimates of the probability of captured via phonon emission. In the two lower 
panels, the positions of the two resonances are indicated. 



